Déplacer des gouttes de liquide peut être très utile, notamment dans les dispositifs micro-fluidiques (lab-on-a-ship) ou afin d’accroitre les performances en transfert thermique (condensation en film ou en goutte…). La modification de la mouillabilité d’un matériau nécessite de mettre en œuvre des traitements physiques (texturation…) et/ou chimiques (greffage, dépôt…) afin de pouvoir ajuster cette propriété de surface, idéalement de superhydrophobe à superhydrophile. Cependant, déplacer une goutte sans action extérieure (gravité, pression…) autre que la mouillabilité nécessite non seulement de pouvoir changer la mouillabilité mais aussi de la faire varier spatialement. La question du profil optimal de mouillabilité de ces matériaux, communément appelés « à gradient de mouillabilité », se pose alors. Optimiser (ici numériquement) une telle propriété nécessite un modèle numérique fiable. Notre choix s’est porté sur l’utilisation de la méthode de Boltzmann sur réseaux (LBM), qui en plus de retrouver l’équation de Navier-Stokes permet de modéliser des problèmes multiphysiques (ici un écoulement diphasique sur une plaque avec prise en compte de la tension de surface et d’une équation d’état représentative du fluide). Après une comparaison des performances (il s’agira ici de maximiser la vitesse de déplacement d’une goutte sur une longueur donnée) de profils de mouillabilité analytiques (quadratique, linéaire…), les premiers résultats d’optimisation (gradient en différences finies + multi-échelle) de cette propriété seront présentés et discutés.